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Baryon-strangeness correlations: a diagnostic of strongly interacting matter 
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The correlation between baryon number and strangeness elucidates the nature of strongly inter- 
acting matter, such as that formed transiently in high-energy nuclear collisions. This diagnostic can 
be extracted theoretically from lattice QCD calculations and experimentally from event-by-event 
fluctuations. The analysis of present lattice results above the critical temperature severely limits 
the presence of qq bound states, thus supporting a picture of independent (quasi) quarks. 

PACS numbers: 25.75.-q, 12.38.Mh, 25.75.Nq, 25.75.Gz 
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The principal goal of high-energy heavy-ion collisions 
is the creation and exploration of a novel state of mat- 
ter in which the quarks and gluons are deconfined over 
distances considerably larger than that of a hadron 
It had originally been assumed that asymptotic freedom 
would cause such matter to behave as a plasma of mass- 
less quarks and gluons interacting with a relatively weak 
screened chromodynamic Coulomb force. While this pic- 
ture is supported qualitatively by the rapid rise in the en- 
tropy density at a temperature of about T c ~ 170 MeV, 
as obtained by lattice QCD calculations, the fact that the 
high-T behavior falls somewhat below that of an ideal 
gas of massless quarks and gluons , indicates that the 
chromodynamic plasma has a more complex structure. 

Indeed, recent results from lattice QCD calculation 
on spectral functions 0, 0, H, [1] suggest the presence 
of bound, color-neutral states well above T c . This has 
led to the suggestion that at moderate temperatures, 
T ~ 1 — 2 T c , the system is composed of medium-modified 
(massive) quarks and gluons together with their (many 
and possibly colored) bound states @, S, 0, E3] ■ 

Furthermore, jet quenching observed at the Relativis- 
tic Heavy-Ion Collider (RHIC) suggests that the matter 
formed is not hadronic [ll|. Moreover, the RHIC col- 
lisions exhibit strong elliptic and radial collective flows 
that are consistent with predictions by ideal fluid dy- 
namics 

IS 

[l3j . If ideal hydrodynamics is indeed the 
correct framework, rapid thermalization must occur and 
this in turn would seem to rule out a weakly interacting 
plasma. Thus, the nature of the matter being created at 
RHIC is still not clarified and the purpose of the present 
Letter is to offer a novel diagnostic tool for eluciding this 
issue by probing the relevant degrees of freedom and their 
correlations. 

In particular we will argue that the correlation between 
the strangeness S and the baryon number B provides a 
useful diagnostic for the presence of strong correlations 
between quarks and anti-quarks. In order to understand 
this, consider first a situation in which the basic degrees 
of freedom are weakly interacting quarks and gluons. 
Then strangeness is carried exclusively by the s and s 
quarks which in turn carry baryon number in strict pro- 
portion to their strangeness, B s = —^S s , thus rendering 



strangeness and baryon number strongly correlated. This 
feature is in stark contrast to a hadron gas in which the 
relation between strangeness and baryon number is less 
intimate. For example, at small baryon chemical poten- 
tial the strangeness is carried primarily by kaons, which 
have no baryon number. 

These elementary considerations lead us to introduce 
the following correlation coefficient, 
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The average (•) is taken over a suitable ensemble of 
events and the last expression uses the fact that (S) van- 
ishes. When the active degrees of freedom arc individual 
quarks, the total strangeness is S = vs — Vs, while the 
baryon number can be expressed as B — ^(U + D — S), 
where U = v u — Vq, is the upness and D = ~ v g is the 
downness. Thus, if the flavors are uncorrelated, we have 
ctbs = — \og and Cbs is unity. 

In a gas of hadron resonances, the total baryon num- 
ber is B = n k B k and its total strangeness is S = 
^2 k n k S k , where the specie k has baryon number Bk and 
strangeness S k . If no prior correlations are present, the 
correlation coefficient may then be expressed in terms of 
the multiplicity variances o\ ~ 



{n k 



Cbs 



Efc (rcfc) B k S k . , 



where the approximate expression holds for Poisson 
statistics. Thus, in the hadronic gas, the numerator 
receives contributions from only (strange) baryons (and 
anti-baryons), while the denominator receives contribu- 
tions also from (strange) mesons, 



Cbs 



(A) + (A) + . . . +3(n-)+3(n+) 
K° + K° + . . . + 9 (Q-) + 9 (n+) 



(3) 



The hadronic freeze-out value value of Cbs is shown in 
Fig. [TJ At the relatively high temperatures relevant at 
RHIC, the strange mesons significantly outnumber the 
strange baryons, so Cbs is smaller than unity. Indeed, 
including hadrons up to f2~, we find Cbs = 0.66 for 
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T = 170 McV and zero chemical potential, ^b = 0. 
As /is is increased, the freeze-out temperature decreases 
[3] and, consequently, strange baryons steadily gain fa- 
vor relative to antikaons, so Cbs increases. As this 
trend continues, only very few antibaryons are present so 
the (positive) strangeness is carried primarily by kaons 
and compensated by A and S hyperons. We then have 
(A + E) ~ (K) and so Cbs ~ §■ This significant depen- 
dence of Cbs orL the hadronic environment is in sharp 
contrast to the simple quark-gluon plasma where the cor- 
relation coefficient remains strictly unity at all tempera- 
tures and chemical potentials. 
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FIG. 1: The correlation coefficient Cbs = o'bs/^s for a 
hadron gas (including all species up to the at freeze- 

out, shown as a function of the baryon chemical potential 
\ib (with the temperature T decreasing from 170 MeV as \ib 
is increased, as found in Ref. [3). Also shown is the corre- 
sponding result for an ideal quark-gluon plasma (dashed line) . 

Whether the chromodynamic system, with respect to 
its baryon-strangeness correlations, does in fact appear as 
an assembly of elementary quasiparticles can be checked 
by lattice calculations. We thus note that in statistical 
equilibrium the correlation coefficient can be expressed 
as a ratio of susceptibilities, Cbs = — ^Xbs/xsSi which 
in turn are second derivatives of the free energy with 
respect to the chemical potential(s), 

1 d 2 F 1 d 2 F 

XBS ~ V d^ B d^s ' XSS ~ ~ V 9/4 ■ ( } 

In terms of the basic flavor densities u, d, s we thus have 
C BS = 1 + Xds + Xus ■ (5) 

Xss 

At fiB = 0, the mixed flavor susceptibilities generally 
tend to be relatively small above T c [l5[. Indeed, using 
the values Xff extracted at T = 1.5 T c by Gavai et al. 
0, we obtain { X us + Xds)/Xss = 0.00(3)/0.53(l) « 1. 
These results were obtained in a quenched approxima- 
tion, but the effect of sea quarks is expected to be 



marginal above T c [lfj. We may thus surmise that 
the lattice system has Cbs ~ 1, suggesting that the 
quark flavors are uncorrelated, as in the ideal quark-gluon 
plasma. (However, the presence of pure gluon clusters 
cannot be ruled out by this diagnostic.) It would obvi- 
ously be of interest to move the lattice calculations be- 
yond the quenched approximation. 

Since the ratio of strange to non-strange susceptibili- 
ties, X ss/(Xuu + Xdd), obtained from lattice QCD agrees 
with that of a hadron gas right at the critical tempera- 
ture T c [l7j|, it would be interesting to see what lattice 
QCD would yield for Cbs at this temperature. Unfortu- 
nately, to our knowledge, no off-diagonal susceptibilities 
involving strange quarks have yet been calculated. How- 
ever, the available results for light quarks [l5| indicate 
that the off-diagonal susceptibilities are smaller than the 
diagonal ones by a factor of twenty at T c . Thus, if flavor 
symmetry holds, the lattice value of Cbs would differ 
from that of the hadron gas already at T c . We also note 
that the vanishing of the off-diagonal susceptibilities, and 
hence the unit value of Cbs, does not conflict with the 
existence of hadron-like resonances that have been iden- 
tified well above T c @, d, @| , since their large masses (of 
more than 2 GeV) make them insignificant near T c . 

It is important to recognize that both the calculated 
equation of state Q and the observed collective flow 
[r3 . [lli | indicate that the system cannot be merely an 
assembly of weakly interacting elementary quarks and 
gluons. However, the various apparently conflicting fea- 
tures might be reconciled if the system were to organize 
itself into an assembly of weakly interacting quasipar- 
ticles, such as the picture emerging from QCD by the 
application of resummation techniques [l8T | . 

Recently there has appeared a model that purports 
to explain both the equation of state as obtained on the 
lattice as well as the large flow observed in heavy-ion col- 
lisions 0, @, Q • The model describes the chromodynamic 
system as a gas of massive quarks, antiquarks, and gluons 
together with a myriad of their bound states generated 
by a screened Coulomb potential. In order to assess the 
consistency of this model with present lattice results, we 
estimate the ratio Cbs m such a scenario. 

We base our estimates on Ref. [1| in which the tempera- 
ture dependence of the screening length and the effective 
masses were obtained by parametrizing lattice results. 
The resulting attraction produces a total of 749 bound 
states, of which only the color-triplet sg and the color- 
singlet qs states (and their conjugates) are of relevance 
here. (The color-hexaplet sg states as well as the diquark 
states are very weakly bound and dissolve entirely at the 
temperature 1.5 T c considered here). There are 4 7r-like 
(spin-singlet) and 12 p-like (spin-triplet) qs states as well 
as 18 sg states (and their conjugates). The abundancies 
of these states are estimated in a grand canonical ensem- 
ble with vanishing chemical potentials. The qs multiplcts 
carry no baryon number and thus contribute only to cr|, 
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hence drive Cbs towards zero, while the more degenerate 
sg multiplet contributes also to gbs an d thus drives Cbs 
towards unity. The resulting value obtained at T = 1.5T C 
is Cbs = 0.62. The fact that it differs significantly from 
the value extracted directly from lattice QCD (see above) 
suggests that the bound-state model, in spite of its in- 
tuitive appeal, cannot account for the BS correlations 
implied by the underlying QCD theory. (No particular 
significance should be attached to the fact that the value 
of Cbs obtained in the bound-state scenario comes out 
as fairly close to the value of a standard hadron gas.) 

In order to estimate what to expect for a RHIC event if 
no QGP is generated, we have extracted Cbs from events 
generated by HIJING [h|] which presents an extension of 
the string fragmentation treatment of hadron production 
in pA and A A collisions. Thus we calculate 



Cbs 
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where and denote the baryon number and 

strangeness observed for a given event n within a central 
rapidity interval, \y\ < 2/ max - The result for Au+Au at 
the top RHIC energy of 100, A GcV is shown in Fig.[2]as a 
function of the rapidity acceptance range y max - At small 
values of y max , only hadrons emerging within a narrow 
slice at mid rapidity arc included and in the limit where 
no event contributes more than one strange hadron the 
resulting value of Cbs equals the ratio of central rapid- 
ity density of baryonic strangeness to the central rapid- 
ity density of all strangeness (The fact that the result, 
Cbs ~ 0.50, lies somewhat below the value shown in 
Fig. [TJ for [j, B = reflects the general underproduction of 
strange baryons by the string mechanism.) As the accep- 
tance is increased, a broader range of chemical potentials 
are sampled and, as a result, the value of Cbs goes up, 
following qualitatively the behavior with fiB exhibited in 
Fig. [TJ Finally, when y max is increased towards the point 
where all baryons are accepted in each event, the baryon 
fluctuation disappears and, consequently, Cbs tends to 
zero (while the strangeness still exhibits a small fluctua- 
tion due to weak decays). 

To illustrate the effect of this canonical suppression 
effect, we include results obtained for e + e~ collisions at 
y/s = 200 GeV with the underlying JETSET code [U, 
for which the net baryon density vanishes at all rapidities. 

As a further reference, the correlation coefficient has 
been exctracted from RHIC events generated by RQMD 
[ID], yielding Cbs ~ 0.55. (This value is also lower 
than that associated with the hadronic freeze-out, re- 
flecting the underestimate of the baryonic strangeness by 
RQMD.) 

While it is relatively straightforward to extract Cbs 
from the various models, as brought out above, it is more 
difficult to determine it experimentally. In principle, Cbs 
can be obtained from event-by-event fluctuation analy- 
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FIG. 2: The correlation coefficient Cbs obtained with HI- 
JING fl9l | from generated samples of 10,000 Au+AU events 
at the top RHIC energy of y/s = 200 A GeV, shown as a 
function of the maximum rapidity accepted, \y\ < y ma x- (The 
points scatter due to the limited event statistics.) Also shown 
is the corresponding results for e + e~ at yfs = 200 GeV cal- 
culated with JETSET [2^]. The arrows point to the values 
obtained for the ideal quark-gluon plasma, the hadronic gas 
at T — 170 MeV, and the bound-state QGP advocated in 
Ref. 1 (at 1.5 T c ). 



sis designed to measure susceptibilities |21| . Since both 



baryon number and strangeness are conserved quantities 
their correlations and fluctuations should be preserved 
during the expansion stage, provided the kinematic ac- 
ceptance is large compared to the typical smearing asso- 
ciated with the hadronization. At the same time it should 
be sufficiently small to render global conservation effects 
negligible. This situation is analogous to the analysis of 
net charge fluctuations made in Rcfs. [22, 23. |24|. 

Relative to that work, the present situation is in fact 
somewhat more favorable, since the principal carriers of 
baryon number and strangeness are heavier than those 
of charge (pions), so one would expect smaller kinematic 
shifts during the hadronization. Furthermore, strong de- 
cays of strange baryons into kaons and nucleons, which 
in principle could wash out the correlation when small 
acceptances are used, are relatively rare 26(. On the 
other hand, an actual experimental investigation should 
correct for the occurrence of weak decays, which remove 
strangeness, and also address the problem of neutron de- 
tection. These issues remain to be studied in detail. 

In this Letter, we have proposed the correlation coef- 
ficient Cbs = —3<tbs/o~s as a useful diagnostic tool for 
elucidating the character of chromodynamic matter. If 
the baryon number and strangeness attributes are car- 
ried by effectively independent quarks and antiquarks, 
then this ratio is unity for any value of the energy and 
baryon- number density. By contrast, it exhibits a sig- 
nificant dependence on the baryon density in a hadron 
gas and is about § at RHIC. We have pointed out that 



4 



the value for Cbs obtained from (quenched) lattice QCD 
calculations are close to unity, while the recently pro- 
posed "bound-state quark-gluon plasma" yields signifi- 
cantly lower values. We also noted that the commonly 
employed event generators HIJING and RQMD yield 
Cbs values at RHIC that are even lower than those asso- 
ciated with the hadronic freeze-out. Thus Cbs appears 
to provide a useful tool for distinguishing between effec- 
tive strangeness-carrying degrees of freedom in various 
models. 

The current experimental activity is concentrated at 
RHIC and we have therefore considered scenarios hav- 
ing vanishing chemical potential, which arc also easier 
to treat with lattice QCD. However, since the hadronic 
Cbs exhibits a marked dependence on /is, it would be 
interesting to explore this observable over a wide range 
of net baryon densities and thus broaden the insight into 
the character of the matter formed. 

Finally, we recall that lattice QCD matter has an 
equation of state that falls below that of a free gas of 
quarks and gluons at high temperature while neverthe- 
less being devoid of flavor correlations (as reflected in the 
vanishing mixed flavor susceptibilities and causing the 
strangeness correlation coefficient Cbs to be unity). We 
wish to suggest that these apparently inconsistent prop- 
erties might be reconciled if strongly interacting matter 
above T c would organize itself into effectively indepen- 
dent but massive quasiparticles. The fact that the lattice 
equation of state can be reproduced in such a quasiparti- 
cle description was already demonstrated in the ud sector 
27j. We emphasize again, however, that the proposed di- 
agnostic is insensitive to the possible existence of strong 
correlations among the gluons (whose presence could lead 
to a hydrodynamic behavior). 

The reduction of the complicated many-body problem 
to the relative simplicity of the quasi-particle picture is 
analogous to the highly successful independent-particle 
description of atomic nuclei and could conceivably 
emerge from the dressing of the elementary quarks and 
gluons by their mutual interaction. In such a picture the 
effective masses arise from a predominantly repulsive 
interaction which in turn would cause a collective 
expansion and thus possibly help explain the observed 
flow in RHIC experiments. 
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